# if (!require('dplyr')) install.packages('dplyr'); library('dplyr')
# if (!require('tidyr')) install.packages('tidyr'); library('tidyr')
if (!require('lubridate')) install.packages('lubridate'); library('lubridate')
if (!require('tidyquant')) install.packages('tidyquant'); library('tidyquant')
if (!require('plotly')) install.packages('plotly'); library('plotly')
Regression Models Testing using Covid Data
Data Ingestion
covid_malaysia_endpoint <- "https://raw.githubusercontent.com/MoH-Malaysia/covid19-public/main/epidemic/cases_malaysia.csv"
covid_malaysia_df <- read.csv(covid_malaysia_endpoint, header=TRUE)
covid_malaysia_df$date <- as.Date(covid_malaysia_df$date, format="%Y-%m-%d")
str(covid_malaysia_df)
## 'data.frame': 705 obs. of 31 variables:
## $ date : Date, format: "2020-01-25" "2020-01-26" ...
## $ cases_new : int 4 0 0 0 3 1 0 0 0 0 ...
## $ cases_import : int 4 0 0 0 3 1 0 0 0 0 ...
## $ cases_recovered : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cases_active : int 4 4 4 4 7 8 8 8 8 8 ...
## $ cases_cluster : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cases_unvax : int 4 0 0 0 3 1 0 0 0 0 ...
## $ cases_pvax : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cases_fvax : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cases_boost : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cases_child : int 0 0 0 0 1 0 0 0 0 0 ...
## $ cases_adolescent : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cases_adult : int 1 0 0 0 2 1 0 0 0 0 ...
## $ cases_elderly : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cases_0_4 : int 0 0 0 0 1 0 0 0 0 0 ...
## $ cases_5_11 : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cases_12_17 : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cases_18_29 : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cases_30_39 : int 0 0 0 0 1 0 0 0 0 0 ...
## $ cases_40_49 : int 1 0 0 0 0 1 0 0 0 0 ...
## $ cases_50_59 : int 0 0 0 0 1 0 0 0 0 0 ...
## $ cases_60_69 : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cases_70_79 : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cases_80 : int 0 0 0 0 0 0 0 0 0 0 ...
## $ cluster_import : int NA NA NA NA NA NA NA NA NA NA ...
## $ cluster_religious : int NA NA NA NA NA NA NA NA NA NA ...
## $ cluster_community : int NA NA NA NA NA NA NA NA NA NA ...
## $ cluster_highRisk : int NA NA NA NA NA NA NA NA NA NA ...
## $ cluster_education : int NA NA NA NA NA NA NA NA NA NA ...
## $ cluster_detentionCentre: int NA NA NA NA NA NA NA NA NA NA ...
## $ cluster_workplace : int NA NA NA NA NA NA NA NA NA NA ...
## [1] 705 31
A brief Look on the graph
fig <- plot_ly(covid_malaysia_df, type = 'scatter', mode = 'lines')%>%
add_trace(x = ~date, y = ~cases_new, name = 'Daily New Cvoid Cases')%>%
layout(showlegend = F)
options(warn = -1)
fig <- fig %>%
layout(
xaxis = list(zerolinecolor = '#ffff',
zerolinewidth = 2,
gridcolor = 'ffff'),
yaxis = list(zerolinecolor = '#ffff',
zerolinewidth = 2,
gridcolor = 'ffff'),
plot_bgcolor='#e5ecf6', width = 1200)
fig
End of Test Document